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Computation  of  Discrete  Slanted  Hole  Film 
Cooling  Flow  Using  the  Navier- Stokes  Equations 

by 


H.  J.  Glbeling,  S.  J.  Shamroth  and  H.  McDonald 
Scientific  Research  Associates,  Inc. 

P.0.  Box  498 
Glastonbury,  CT  06033 


ABSTRACT 


An  analysis  and  computational  procedure  have  been  developed  for  predicting 
flow  and  heat  transfer  which  results  from  coolant  injection  through  a  single 
row  of  round  holes  oriented  at  an  angle  to  a  flat  surface  with  the  injection 
and  free  stream  velocity  vectors  coplanar.  This  method  solves  the  compressible 
Navier-Stokes  equations  and  utilizes  "zone  embedding",  surface-oriented  co¬ 
ordinates,  interactive  boundary  conditions,  and  an  efficient  split  LBI  scheme. 

The  approach  treats  the  near-hole  flow  region  where  the  film  cooling  flow  is 
initially  established.  Prior  studies  considered  only  laminar  flow  in 
order  to  simplify  development  of  the  computational  procedure.  Under  the  present 
effort,  several  turbulence  models  suitable  for  the  discrete  hole  film  cooling 
problem  have  been  Investigated  by  predicting  a  number  of  test  cases  from  the 
1980-81  AFOSR-HTTM- S tanf o rd  Conference  on  Complex  Turbulent  Flows.  The  cal¬ 
culation  of  a  discrete  hole  film  cooling  case  for  an  injection  angle  of  35  degrees 
has  been  initiated  using  a  mixing  length  turbulence  model.  These  results  will 
be  compared  with  the  experimental  data  of  Kadotani  and  Goldstein,  and  recal¬ 
culation  with  either  a  turbulence  kinetic  energy/algebraic  length  scale  or  a 
turbulence  kinetic  energy/dissipation  rate  model  would  be  carried  out  under 
the  subsequent  effort,  if  warranted  by  the  data  comparison. 
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INTRODUCTION 

To  achieve  higher  turbine  efficiencies,  designers  have  been  forced  to 
deal  with  high  turbine  inlet  tenperatures  which  may  exceed  structural  limita¬ 
tions  of  available  materials  used  for  turbine  blades  unless  some  method  of 
cooling  these  blades  is  Introduced.  Various  cooling  schemes  have  been  devised 
to  ensure  that  the  structural  integrity  of  the  blades  remain  Intact  when 
.  exposed  to  a  high  temperature  environment  for  prolonged  periods  of  time.  One 
of  the  more  promising  schemes  is  discrete  hole-  film  cooling,  wherein  cooling 
air  is  injected  through  either  a  single  row  or  multiple  staggered  rows  of  holes, 
to  provide  a  protective  layer  of  cool  air  between  the  blade  surface  and  hot 
mainstream  gases.  Discrete  hole  film  cooling  is  better  suited  to  the  high 
temperatures  of  advanced  high  performance  engines  than  the  current  convectively 
cooled  blades  found  in  commercial  engines.  A  large  number  of  parameters  can  be 
varied  to  obtain  a  configuration  which  optimizes  heat  transfer  and  aerodynamic 
characteristics,  including  hole  shapes  and  patterns,  injection  angles  and 
rates,  surface  curvature,  coolant  temperatures,  and  mainstream  boundary  layer 
characteristics.  Since  experimental  determination  of  optimal  configurations 
is  costly,  a  computational  procedure  which  could  be  used  to  screen  alternative 
configurations  without  experimental  testing  would  be  of  considerable  value. 

An  analysis  aid  computational  procedure  were  developed  previously  for  predicting 
flow  and  heat  transfer  which  results  from  coolant  injection  through  a  single  row 
of  round  holes  oriented  at  a  non-normal  angle  to  a  flat  surface.  This  work  was 
an  extension  of  a  previous  effort  which  considered  only  normal  injection 
through  a  single  row  of  round  holes.  Both  of  these  prior  studies  considered 
only  laminar  flow  in  order  to  simplify  development  of  the  computational  procedure. 
The  present  effort  involves  investigation  of  turbulence  models  suitable  for 
the  discrete  hole  film  cooling  problem.  Therefore,  the  model  chosen  must  be 
capable  of  treating  flows  with  large  pressure  gradients  and  separated  regions. 
Several  turbulence  models  were  investigated  under  the  present  study  by 
predicting  a  number  of  test  cases  from  the  1980-81  AFOSR-HTTM-Stanford 
Conference  on  Complex  Turbulent  Flows.  The  models  considered  were: 

*Two  of  the  cases  reported  on  herein  were  partially  funded  by  AFOSR 
and  the  U.S.  Army  Research  Office. 
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(1)  algebraic  mixing  length  (1),  (2)  turbulence  kinetic  energy  and  algebraic 
length  scale  (k-A),  and  (3)  turbulence  kinetic  energy  and  dissipation  rate  (k-e). 
All  of  these  models  utilize  an  isotropic  eddy  viscosity  formulation  to  represent 
the  Reynolds  stress  terms.  This  study  will  ultimately  provide  an  Improved 
understanding  of  the  turbulent  flow  film  cooling  process  *  particularly  in  the 
near-hole  region,  and  represents  the  third  step  toward  development  of  a  computer 
program  capable  of  predicting  heat  transfer  and  aerodynamic  loss  levels 
associated  with  discrete  hole  film  cooling  on  gas  turbine  blades. 


PREVIOUS  WORK 


Previous  experimental  work  on  film  cooling  flow  has  consisted  mainly  of 
flow  visualization  studies  and  measurements  of  film  cooling  effectiveness  down¬ 
stream  of  the  coolant  injection  hole  for  a  variety  of  flow  conditions.  For 
example,  Colladay  and  Russell  [1]  performed  a  flow  visualization  study  of  dis¬ 
crete  hole  film  cooling  in  three  different  hole  configurations,  to  obtain  a 
1  better  understanding  of  flow  behavior  associated  with  coolant  injection.  They 
considered  injection  from  discrete  holes  in  a  three-row  staggered  array  with 
five-diameter  spacing  and  three-hole  angles:  (1)  normal  to  the  surface, 

(2)  slanted  30  degrees  to  the  surface  but  aligned  (coplanar)  with  the  free 
stream,  and  (3)  slanted  30  degrees  to  the  surface  and  45  degrees  laterally  to 
the  free  stream.  For  flow  conditions  typical  of  gas  turbine  applications, 

Colladay  and  Russell  observed  that  normal  injection  is  subject  to  flow  separa¬ 
tion  even  at  low  injection  rates,  that  slanted  injection  works  well  except  at 
high  Injection  rates,  and  that  compound  slanted-yawed  Injection  is  less  susepti- 
ble  to  separation  but  tends  to  promote  entrainment  of  free  stream  fluid  (which 
augments  heat  transfer)  and  to  increase  secondary  flow  losses.  These  experimental 
observations  emphasize  the  importance  of  Interaction  between  the  Injected  fluid 
and  the  free  stream. 

Kadotani  and  Goldstein  [2]  have  performed  experimental  measurements  for 
discrete  hole  film  cooling  through  a  row  of  holes  having  three-diameter  lateral 
spacing.  Inclined  at  35  degrees  to  the  injection  surface,  and  coplanar  with  the 
free  stream.  They  considered  a  variety  of  flow  conditions  and  varied  boundary 
layer  thickness  and  injection  rates,  as  well  as  Reynolds  number  and  turbulence 
properties.  They  measured  lateral  and  streamwise  distributions  of  film  cooling 
effectiveness  downstream  of  injection,  and  also  measured  mean  velocity  and 
temperature  in  a  transverse  plane  normal  to  the  free  stream  and  located  at  the 
downstream  edge  of  the  injection  holes.  Kadotani  and  Goldstein  deduced  from 
their  measurements  that  reversed  flow  is  present  near  the  injection  hole  even 
for  injection  rates  as  low  as  0.35,  and  that  the  size  and  strength  of  the 
reversed  flow  Increases  with  increasing  boundary  layer  thickness.  They  also 
found  that  the  approaching  boundary  layer  thickness  significantly  affects  the 
flow  behavior  and  lateral  variation  of  film  cooling  effectiveness,  particu¬ 
larly  near  the  injection  hole. 


In  a  prior  study  of  injection  through  a  single  array  of  round  holes 
oriented  normal  to  a  flat  surface,  Kreskovsky,  Briley  and  McDonald  [3]  solved 
the  compressible  Navler-Stokes  equations  utilizing  surface-oriented  coordinates 
and  an  efficient,  consistently-split  LBI  scheme.  Subsequently,  Gibeling, 
Kreskovsky,  Briley  and  McDonald  [A]  extended  the  work  of  [3]  to  treat  the  coolant 
injection  through  a  single  row  of  round  holes  oriented  at  a  non-normal  angle  to 
a  flat  surface.  In  this  analysis  the  governing  equations  in  Cartesian  coordinates 
‘  were  cast  into  strong  conservation  form  using  a  Jacobian  transformation  to 
obtain  the  governing  equations  in  a  general  nonorthogonal  coordinate  system. 

Both  previous  studies  [3  and  A]  focused  on  the  near-hole  flow  region  where  the 
film  cooling  flow  is  initially  established  and  where  local  hot  spots  may  occur. 
Since  the  flow  field  surrounding  injection  may  involve  separation  and  reversed 
flow  relative  to  the  free  stream,  and  since  directional  and/or  other  properties 
of  the  flow  which  might  permit  simplifying  assumptions  are  not  clearly  evident, 
these  analyses  are  based  on  numerical  solution  of  the  three-dimensional 
compressible  Navier-Stokes  equations.  In  selecting  the  computational  domain, 
a  "zone  embedding"  approach  is  adopted  by  [3, A]  whereby  attention  is  focused  on 
a  subregion  of  the  overall  flow  field  in  the  immediate  vicinity  of  the  discrete 
hole  coolant  injection.  At  the  curved  boundaries  located  within  the  free  stream 
region,  interactive  boundary  conditions  which  permit  inflow  and  outflow  of  shear 
layers  and  the  compressible  inviscid  free  stream  are  derived  from  an  assumed 
flow  structure.  Since  details  of  the  coolant  velocity  distribution  at  the 
hole  exit  are  not  known  a  priori  and  presumably  will  depend  on  an  interaction 
between  the  coolant  and  the  free  stream  flow,  the  computational  domain  is 
chosen  to  include  the  flow  region  within  the  coolant  hole  as  well  as  that 
exterior  to  the  hole.  Including  the  flow  region  within  the  coolant  hole 
permits  the  interaction  between  the  coolant  and  free  stream  flows  and  its 
important  influence  near  the  hole  to  be  determined  as  part  of  the  final 
solution,  without  the  need  for  simplifying  assumptions.  Sample  laminar  cal¬ 
culations  yielded  predictions  of  the  interaction  between  the  injectant  and 
main  stream  flow  which  were  in  qualitative  agreement  with  experimental 
observations. 

The  only  other  previous  computational  study  known  to  the  authors  which 
treated  the  discrete  hole  film  cooling  problem  using  a  three-dimensional  cal¬ 
culation  procedure  is  that  of  Bergeles,  Gosman  and  Launder  [5],  who  applied 
the  approximate  "partially  parabolic"  calculation  procedure  of  Pratap  & 

Spalding  [6]  for  the  case  of  laminar  flow.  This  calculation  procedure  neglects 
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streamwise  diffusion  and  employs  iterated  forward  marching  solution  of  three 
approximate  momentum  equations.  The  procedure  begins  with  a  "guessed"  pressure 
field  and  performs  Iterated  forward  marching  sweeps  of  the  three-dimensional 
flow  field,  solving  the  approximate  momentum  equations  and  utilizing  various 
strategies  to  modify  or  correct  the  pressure  field,  so  as  to  improve  the 
continuity  balance.  Since  the  calculation  procedure  is  based  on  forward  marching 
the  flow  region  within  the  coolant  injection  hole  is  excluded  from  the  computa¬ 
tional  domain,  and  instead,  the  coolant  velocity  distribution  is  specified  at 
the  hole  exit  as  a  two-dimensional  uniform  stream  at  a  prescribed  Injection 
angle.  Bergeles,  Gosman  and  Launder  made  flow  calculations  for  several  flow 
conditions,  including  different  injection  angles. 


TURBULENCE  MODELS 


This  section  presents  the  turbulence  model  equations  in  the  general 
form  utilized  in  the  test  cases  considered  herein.  Details  concerning 
algebraic  length  scale  specification  are  given  in  the  sections  describing 
each  test  case,  since  the  length  scale  formulation  must  be  adapted  to  the 
particular  flow  being  computed.  The  latter  requirement  is  a  major  disadvantage 
of  models  using  a  specified  length  scale;  therefore,  consideration  is  given 
to  the  two-equation  (k-e)  model  [7,8]  with  appropriate  low  turbulence  Reynolds 
number  corrections  [8]  to  permit  sublayer  resolution. 

Turbulence  Kinetic  Energy  Equation 

Several  cases  presented  in  this  report  utilize  a  turbulence  kinetic 
energy  equation  and  an  algebraic  length  scale  equation.  The  turbulence  kinetic 
equation  as  given  by  [7,8]  is 

^■+V  (/,lTk)  =  V-[(M  +  ^)Vk]  +/i7(2&:l5) -P€ -2Ai(Vk,,z)  (Vk,,JJ)  (1) 


where  k  is  the  turbulence  kinetic  energy. 


Turbulence  Energy  Dissipation  Equation 


(2) 


The  equation  for  the  dissipation  rate  of  turbulence  energy,  e,  as 
given  by  [7,8]  is 


V(pvTt)  =  V-[  (^  + 


(3) 


-  *  /  a2ir  \  /  a2»r  \ 


where  the  deformation  tensor  D  is 


D=  ^-[(Vu)+(7Tr)T]  (4) 

The  turbulent  flow  stress  tensor  (Reynolds  stress)  which  appears  in  the 
momentum  conservation  equations  [9,10]  is  modeled  using  an  isotropic  eddy 
viscosity  formulation,  i.e., 

irT  -  -  pu‘  u*  =  2 "5  -  fJ-jV *"u  (5) 

s< 

The  turbulent  viscosity  yt  appearing  in  Eqs.  (1,3,5)  must  be  determined  using 
a  suitable  model,  and  the  parameters  o^,  o£,  and  C2  must  be  specified. 

Turbulence  Model  Relations 

The  turbulent  viscosity  introduced  above  is  obtained  from  the  Prandtl— 
Kolmogorov  relation  viz. 

pkz 

H  S  CM  ~T~  (6) 


When  the  dissipation  rate  equation  (3)  is  not  solved ,  e  is  determined  from 
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c  3/4  k3/2 

*  T 


(7) 


where  the  turbulence  length  scale,  1,  must  be  specified  consistent  with  the 
expected  turbulence  structure  of  the  flow. 

The  constraints  ok,  a£  and  are  taken  as  [7], 


°k  *  ' 0 

crf  r  1.3 

C,  =  1.44 

for  fully  developed  turbulet  flow  C  -0.09  and  C2 
flows  Jones  and  Launder  [8]  recommend. 


(8) 

1.92;  for  relaminarizing 


=  0.09  exp [-2. 5 (1.0+  RT/50)j 
C2  =  2.0 [  1.0  -  0.3  exp(  — RT2)]  (9) 


However,  these  expressions  have  not  given  successful  predictions  of  flows 
undergoing  forward  transition  (e.g.,  [11]).  An  alternate  transition  model 
was  developed  by  Shamroth  and  Gibellng  [10]  and  was  successfully  applied  to 
the  time-dependent  airfoil  flow  field  problem.  This  method  used  the  turbulence 
kinetic  energy  equation  excluding  the  low  Reynolds  number  correction  term 
(last  term  in  Eq.  (1)  ),  in  conjunction  with  an  algebraic  length  scale  relation. 


The  analysis  of  [10]  follows  the  Integral  turbulence  energy  procedure  of  [12-14] 
by  introducing  a  turbulence  function  a^, 

°l  *  Cftl'Z  /Z  (10) 


and  a^  is  taken  as 


function  of  a  turbulence  Reynolds  number,  R_ ,  of  the  form 


where 


+  6.66  o0 


f(Rr) 

100 


(11) 


o0  -  .0115 

f(Rr)  ■  100.  Rt022  Rr  <  I 
f(RT)  •  68  I  Rt  +  614.3  Rt  >40 


(12) 


and  a  cubic  curve  fit  was  used  for  values  of  between  1  and  40.  Since 
[12-14]  used  an  integral  procedure  an  average  value  of  RT  was  defined.  In 
[10],  a  local  value  of  Rt  was  defined  as 


AS 


(13) 


Then  a^  was  evaluated  via  Eq.  (11)  and  follows  from  Eq.  (10).  The 
transition  model  of  [10]  was  utilized  in  several  of  the  cases  presented  herein; 
each  with  an  appropriately  specified  algebraic  length  scale  relation. 


Turbulence  Length  Scale 


An  algebraic  length  scale  relation  is  required  in  the  turbulence  model 
utilizing  only  the  turbulence  kinetic  energy  equation.  The  near  wall  region 
nrfT-tng  length  is  obtained  from  the  McDonald-Camarata  model  [15]  with  Van  Driest 
damping  [16] , 


l  -  lm ton  h 


*00 


(14) 


where  k  is  the  von  Karman  constant  and  A 


is  the  van  Driest  damping  coefficient. 


K  =  0.4 
A*  -  26.0 


(15) 


and  l 


096,  with  the  boundary  layer  thickness  6  to  be  specified. 


The  nondlmenslonal  distance  y  is  defined  as 


end  the  friction  velocity  U  in  the  present  analysis  is  taken  as 


(16) 


(17) 


where  the  local  shear  stress  is  obtained  from 

=  /t.ff(2D:D),,Z  (18) 

Note  that  for  small  y  the  tanh  function  in  Eq.  (14)  reduces  to  icy  while  for 
large  y  it  approaches  lB. 

A  more  detailed  description  of  the  actual  turbulence  models  used  for  the 
computation  and  their  effect  on  the  computation  is  given  in  the  results 


section. 


RESULTS 


As  part  of  the  present  effort,  several  test  cases  were  computed  from  the 
1980-81  AFOSR-HTTM- Stanford  Conference  on  Complex  Turbulent  Flows.  Case  number 
0423,  the  primary  case  computed  under  this  effort,  considers  the  incompressible 
flow  over  a  backward-facing  step  with  a  10°  turning  downstream  of  the  step. 

This  was  a  predictive  case  for  which  the  experimental  results  were  not  provided 
until  after  the  conference.  Case  numbers  8601  and  8641  were  funded  primarily 
by  the  U.S.  Army  Research  Office  with  computer  resources  funded  by  AFOSR.  Case 
8601  considers  the  compressible  flow  in  a  tube  containing  a  normal  shock  wave. 
The  data  for  this  case  was  taken  by  Mateer,  Brosh  and  Viegas  [17].  Case  8641 
considers  a  planar  free-shear  layer  generated  by  a  backward  facing  step, 
reattacking  on  a  20°  ramp  at  supersonic  speed.  The  data  for  this  case  was 
taken  by  Settles,  Baca,  Williams  and  Bogdonoff  [18]. 

Stanford  Case  0423 

A  schematic  of  the  geometry  for  this  case  is  presented  in  Fig.  1. 

A  coordinate  system  (Fig.  2)  for  this  geometry  was  developed  under  the  present 

effort  beginning  with  Cartesian  coordinates  upstream  of  the  step,  followed  by 
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a  transition  region  in  which  straight  transverse  coordinate  lines  (y  -direction) 
are  rotated  through  10  degrees,  and  ending  with  another  Cartesian  region  at  the 
final  10-degree  angle. 

Two  turbulence  models  were  Investigated  using  this  test  case.  First,  the 
two-equation  (k-e)  model  with  the  Jones-Launder  C^(ut>  relationship  was 
utilized.  Preliminary  computations  were  carried  out  for  the  channel  upstream 
of  the  step  in  Fig.  1,  in  order  both  to  verify  the  k-e  formulation  in  the 
computer  code  and  to  obtain  consistent  upstream  flow  conditions  for  the 
backstep  problem.  In  the  straight  channel  flow  application  the  k-e  model, 

Eqs.  (1-6,  8-9)  gave  quantitatively  accurate  predictions  for  boundary  layer 
velocity  profiles,  skin  friction,  and  flow  development.  However,  when  the 
model  was  applied  to  the  backstep  configuration  numerical  difficulties  in 
solving  the  equations  were  encountered  and  converged  solutions  could  not  be 
obtained.  These  problems  were  found  to  be  directly  related  to  the  low 
turbulence  Reynolds  number  correction  terms  in  Eq.  (1,3).  These  terms  are 
ad  hoc  differential  terms  which  were  devised  by  [8]  for  an  attached 
relamlnarizing  boundary  layer.  In  the  separated  shear  layer  which  occurs 
downstream  of  the  step  the  shearing  tens  (l.e.  2D: B)  Increase  substantially 


with  a  concurrent  Increase  in  both  k  and  e  in  that  region.  The  large 
gradients  in  k,  c  and  the  velocities  lead  to  large  low  Reynolds  number 
correction  terms.  This  results  in  nonphysical  behavior  such  as  excessively 
large  turbulent  viscosities  in  certain  areas  of  the  flow,  and  negative  k  and  e 
in  some  areas.  Numerous  changes  in  the  numerical  treatment  of  these  terms 
gave  little  improvement  in  the  predictions,  and  several  ad  hoc  modifications 
to  the  low  Reynolds  number  correction  terms  did  not  Improve  the  results. 

Therefore,  the  k-e  model  was  abandoned  for  this  case,  and  a  turbulence 
kinetic  energy  model  with  an  algebraic  length  scale  relation,  Eqs.  (1,4-6, 7-8, 
10-18),  was  implemented.  The  mixing  length  distribution  was  constructed  by 
using  Eqs.  (14-18)  with  the  known  boundary  layer  thicknesses  upstream  of  the 
step,  and  an  estimated  thickness  beyond  the  flow  reattachment  point  on  the 
lower  wall  in  Fig.  1.  In  the  region  just  downstream  of  the  step,  the  length 
scale  was  varied  smoothly  from  its  upstream  value  to  its  downstream  value. 

The  computed  flow  field  for  this  case  exhibited  a  cyclic  behavior  In  time  and 
an  absolute  convergence  to  a  steady  solution  could  not  be  obtained.  It  is 
not  known  if  this  is  physical  behavior,  or  if  it  is  caused  by  the  turbulence 
model  or  the  computational  grid  employed.  Furthermore,  the  predicted  re¬ 
attachment  point  was  about  seventeen  step  heights  downstream  of  the  step 
which  is  considerably  more  than  the  experimentally  observed  .value  of  9.4 
step  heights  [19].  A  velocity  vector  plot  for  the  predicted  flow  is  shown 
in  Fig.  3  at  the  last  time  step  in  the  calculation.  These  results  show  an 
apparent  anomalous  behavior  in  the  separated  region,  where  the  flow  appears 
to  be  reattaching  but  then  the  separation  region  enlarges  and  continues  much 
further  downstream  before  actually  reattaching.  The  reasons  for  this  behavior 
are  still  under  investigation;  however,  it  may  be  related  to  the  computational 
grid  variation  in  that  region  if  it  is  a  nonphysical  phenomena. 

Stanford  Case  8601 

Case  8601  [17],  which  considers  the  compressible  flow  in  a  tube  containing 
a  shock  wave,  was  computed  by  Roscoe  and  McDonald  [20] .  Although  the  geometry 
for  this  case  is  simple,  the  fluid  mechanics  is  complicated  because  the 
stationary  shock  position  in  the  tube  is  determined  only  by  the  boundary  layer 
growth  and  the  shock  wave-boundary  layer  Interaction.  Thus,  the  physics 
of  the  near-wall  turbulence  phenomena  must  be  modeled  sufficiently  well  to 


obtain  Meaningful  results.  This  calculation  was  initiated  with  the  k-1 
turbulence  model  described  above  with  zero  k  Imposed  as  the  wall  boundary 
condition.  The  length  scale  specification  is  straightforward  as  long  as  the 
boundary  layer  remains  attached  in  the  interaction  region.  After  the  predicted 
flow  field  was  well  established  in  the  k-1  model  calculation,  the  K-e  model 
described  above  was  introduced  and  the  calculation  was  run  to  convergence. 

The  boundary  conditions  for  k  and  e  were  both  zero  function  values  at  the 
solid  wall.  Other  boundary  conditions  used  in  the  calculation  are  detailed  in 
[20].  Also,  a  time-dependent  adaptive  mesh  in  the  axial  direction  was 
implemented  by  [20]  in  order  to  properly  resolve  the  shock  at  all  times 
during  the  computation.  The  final  computational  mesh  for  this  case  is  shown 
in  Fig.  4  with  the  mesh  lines  near  the  tube  wall  omitted  for  clarity.  There 
were  a  total  of  41  mesh  points  in  the  radial  direction  with  a  high  concentration 
near  the  wall  for  turbulent  boundary  layer  resolution,  and  31  mesh  points 
in  the  axial  direction  with  a  concentration  centered  around  the  shock  location. 

Pressure  contours  in  the  vicinity  of  the  shock  are  shown  in  Fig.  5. 

As  expected,  the  shock  is  sharp  in  the  core  region  and  diffuses  as  the  wall 
is  approached  due  to  the  boundary  layer  interaction.  In  Fig.  6  the  predicted 
wall  pressure  distribution  is  compared  with  the  experimental  results,  and  the 
agreement  is  seen  to  be  quite  good.  A  velocity  vector  plot  from  a  radial 
location  of  0.5R  to  the  tube  wall  is  shown  in  Fig.  7.  This  figure  clearly 
displays  the  shock-boundary  layer  interaction  which  results  in  boundary  layer 
thickening  and  deflection  of  the  outer  flow,  followed  by  recovery  of  the 
boundary  layer  downstream  of  the  interaction  region.  The  measured  and  predicted 
velocity  profiles  are  compared  in  Fig.  8,  where  the  agreement  is  quite  good. 

The  largest  difference  occurs  at  a  streamwlse  location  z/$u  “  8,  which  is 
in  the  recovery  region  fairly  far  downstream  of  the  interaction  region. 

Finally,  the  Reynolds  stress  profiles  are  compared  in  Fig.  9,  where  larger 
discrepancies  between  the  experiment  and  predictions  are  observed.  However, 
it  is  noted  that  data  exhibits  nonphysical  behavior  near  the  interaction  region, 
probably  due  to  probe  interference  or  measurement  difficulties  in  the  transonic 
flow  regime  with  a  hot  wire  probe.  Physically  one  expects  the  Reynolds  stress 
to  peak  in  the  shear  layer,  and  the  location  of  that  peak  is  not  expected  to 
shift  nearly  as  much  as  indicated  by  the  data.  Although  the  two-equation 
turbulence  aodel  employed  gives  good  mean  flow  predictions  for  this  case,  the 
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apparent  problems  with  the  Reynolds  stress  measurements  do  not  permit  any 
further  evaluation  of  the  model  to  be  made. 

Stanford  Case  8641 

Case  8641  [18]  considers  a  planar  free-shear  layer  generated  by  a  backward 
facing  step  with  reattachment  on  a  20°  ramp  at  supersonic  speed.  The  predic¬ 
tions  presented  here  are  taken  from  Weinberg,  et  al  [21].  Settles,  et  al  [18] 
carefully  adjusted  the  experiment  so  that  the  free  shear  layer  left  the  back 
step  parallel  to  the  wall  without  the  expansion  waves  and  subsequent  recom- 
presslon  shock  which  would  occur  for  a  simple  back  step  flow.  As  a  result, 
the  geometry  (Fig.  10)  which  must  be  considered  is  more  complicated.  The 
computational  domain  developed  by  [21]  for  this  case  is  shown  in  Fig.  11.  The 
corner  at  the  step  was  kept  sharp  since  it  has  a  significant  influence  on  the 
flow  field  development,  while  the  corner  at  the  cavity-ramp  juncture  was 
rounded,  since  it  was  not  expected  to  influence  the  reattachment  appreciably. 

The  flow  conditions  for  this  case  are  M  «  2.92,  Re  -  6. 7x10^ /meter, 

40  x 

and  Tw  *  265  K.  The  boundary  conditions  for  this  problem  are  discussed  by  [21). 
The  turbulence  model  used  was  an  algebraic  mixing  length  model  with  the  free 
stream  length  scale  lw  linearly  increased  in  the  wake  region  downstream  of  the 
step  until  it  reached  a  value  of  ■  0.096aVg,  where  6^^  is  the  average 
boundary  layer  thickness  on  the  ramp.  This  resulted  in  a  value  la,  m  0.039  for 
the  present  case. 

Predicted  pressure  contours  are  shown  in  Fig.  12.  The  compression  waves 
generated  by  the  reattachment  coalesce  into  a  shock  as  expected.  The  calculated 
velocity  vectors  presented  in  Fig.  13  show  the  shear  layer  parallel  to  the 
cavity  wall,  the  presence  of  three  vortices  in  the  cavity,  and  the  reattachment 
shock.  Velocity  profiles  in  the  shear  layer  and  on  the  ramp  are  compared  with 
experimental  data  in  Figs.  14  and  15.  The  agreement  between  the  calculations 
and  experiment  is  fairly  good.  The  discrepancies  which  are  observed  are  due 
in  part  to  the  mixing  length  turbulence  model  employed,  the  close  proximity 
of  the  computation  exit  plane  to  the  interaction  region,  and  the  lack  of  mesh 
refinement  near  the  ramp  wall  [21]. 


Turbulent  Film  Cooling  Flow 


Under  the  present  effort  a  comparison  with  the  turbulent  film  cooling 
experiment  of  Kadotani  and  Goldstein  [2]  has  been  initiated.  Initially  a 
mixing  length  turbulence  model  will  be  utilized.  As  part  of  the  follow-on 
effort  this  case  will  be  considered  in  more  detail  using  both  a  turbulence 
kinetic  energy-specified  length  scale  model  (k-1) ,  and  a  two-equation  model 
(k-e).  The  objective  of  this  study  is  determination  of  an  adequate  turbulence 
model  to  permit  accurate  calculation  of  the  film  cooling  flow  field  and  a 
better  understanding  of  the  factors  affecting  the  film  cooling  process. 
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Coordinate  system  for  Stanford  Case  0423  geometry 
Certain  coordinate  lines  omitted  for  clarity. 
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Computational  mesh  for  Stanford  Case  8601:  41  radial  lines  by 

31  axial  lines.  6  is  the  upstream  boundary  layer  thickness. 
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Figure  6  -  Wall  pressure  distribution,  Stanfi 
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Figure  9  -  Reynolds  stress  profiles  at  selected  axial  locations,  Stanford  Case  8601. 
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Figure  14  -  Velocity  Profiles  Through  Shear  Layer. 


